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SUMMARY 


A  method  for  the  prediction  of  the  response  of  a  wing  to  non¬ 
stationary  buffet  loads  is  presented.  The  time  history  of  the  applied  load  is 
segmented  into  a  number  of  time  intervals.  In  each  time  segment,  the  non¬ 
stationary  load  is  represented  by  the  product  of  a  deterministic  shaping 
function  and  a  statistically  stationary  random  function.  An  approximate 
modelling  of  the  load  on  the  wing  is  given.  The  wing  is  divided  into  panels 
or  elements,  and  the  load  is  computed  from  measured  or  estimated  pressure 
fluctuations  at  the  centre  of  each  panel.  A  series  representation,  with  terms 
of  the  correlated  noise  type,  is  used  to  curve  fit  the  experimentally  deter¬ 
mined  complex  buffet  pressure  power  spectral  densities.  Using  the  correlated 
noise  form  of  power  spectral  density  for  the  random  part  of  the  applied 
load,  analytic  expressions  are  derived  for  the  mean  square  displacement  and 
acceleration  response  of  the  wing.  An  illustration  using  data  available  for 
the  F-4E  aircraft  is  included. 


RESUME 


La  presente  communication  porte  sur  une  methode  de  prevision 
du  comportement  d’une  aile  sous  des  charges  de  buffeting  instationnaires. 
L ’evolution  de  la  charge  appliquee  est  divisee  en  intervalles  de  temps.  Dans 
chaque  intervalle,  la  charge  non  fixe  est  representee  par  le  produit  d’une 
fonction  de  “forme”  determinate  et  d’une  fonction  aleatoire  statistiquement 
stationnaire.  Un  modele  approche  de  la  charge  sur  l’aile  est  donne.  L’aile 
est  divisee  en  panneaux  ou  elements,  et  la  charge  est  calculee  a  partir  des 
fluctuations  de  pression  mesurees  ou  estimees  au  centre  de  chaque  panneau. 
Un  developpement  en  serie,  comportant  des  termes  du  type  bruit  correle, 
sert  a  ajuster  a  une  courbe  les  densites  spectrales  complexes  de  pression  de 
buffeting  qui  sont  determinees  experimentalement.  Des  expressions  analy- 
tiques  du  carre  moyen  de  la  reponse  de  1’aile  en  deplacement  et  en 
acceleration  sont  derivees  de  la  fonction  de  type  bruit  correle  de  la  densite 
spectrale  pour  la  composante  aleatoire  de  la  charge  appliquee.  Un  exemple 
s’appuyant  sur  les  donnees  existantes  sur  1 ’avion  F-4E  est  donne. 
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A  METHOD  FOR  THE  PREDICTION  OF 
WING  RESPONSE  TO  NON  -STATION  ARY  BUFFET  LOADS 

1.0  INTRODUCTION 

Buffeting  is  the  aeroelastic  response  of  aircraft  structures  to  aerodynamic  excitation  arising 
from  random  loading  due  to  flow  separations  on  the  wing.  It  is  almost  always  encountered  when  the 
aircraft  approaches  the  limiting  useable  lift  at  high  speeds.  The  maneuvering  capability  of  aircrafts  in 
the  transonic  speed  range  is  thus  usually  limited  by  buffet  or  buffet-related  unsteady  phenomena  which 
induce  the  pilot  to  restrict  the  maneuver.  Methods  for  predicting  the  buffet  intensity  as  the  aircraft 
penetrates  into  the  buffet  regime  are  extremely  useful  and  much  needed  in  aircraft  design. 

The  random  nature  of  the  loading  on  the  wing  due  to  flow  separations  requires  statistical 
theory  in  predicting  the  dynamic  response  of  the  wing  during  buffeting.  Considering  the  exciting  force 
to  be  statistically  stationary,  Liepman  (Ref.  1)  examined  the  problem  of  the  lift  force  exerted  on  a  two- 
dimensional  thin  airfoil  moving  in  turbulent  air.  Later,  he  extended  the  method  to  wings  of  finite  span 
(Ref.  2).  The  analysis  was  generalized  by  Ribner  (Ref.  3)  using  a  model  of  turbulence  represented  by  the 
superposition  of  plane  sinusoidal  shear  waves  of  all  orientations  and  wavelengths.  The  correlations  be¬ 
tween  flight  and  wind  tunnel  tests  with  predictions  based  on  statistical  approach  had  been  reported  by 
various  investigators  (e.g..  Refs.  4-8).  Amongst  the  more  recent  studies  on  this  subject  are  those  of 
Mullans  and  Lemley  (Ref.  9),  Hwang  and  Pi  (Refs.  10, 11).  Cunningham  et  al.  (Ref.  12),  Jones  (Ref.  13), 
and  Butler  and  Spavins  (Ref.  14).  Because  of  the  complexities  of  the  buffet  phenomenon,  the  above 
methods  incorporate  numerous  simplifying  assumptions  and  hence,  can  only  be  considered  to  be  approx¬ 
imate.  However,  they  are  invaluable  in  providing  buffet  information  during  the  early  stages  of  aircraft 
design. 


The  assumption  of  statistical  stationarity  used  in  the  above  prediction  methods  can  be  justified 
as  long  as  the  flight  conditions  are  unchanged  within  a  reasonably  long  analysis  time  when  measurements 
are  made  of  either  the  buffet  load  or  the  wing  response.  In  wind  tunnel  testing,  statistical  stationarity 
requires  long  running  time  compared  to  the  natural  period  of  the  mode  of  vibration  under  investigation. 
When  these  conditions  are  not  met,  such  as  in  rapid  maneuvers  or  in  short  duration  wind  tunnel  testing, 
a  non-stationary  analysis  of  the  wing  response  has  to  be  considered. 

An  approximate  solution  of  the  transient  response  of  a  wing  to  a  buffet  load  having  a  parabolic 
variation  of  its  rms  value  with  time  was  given  by  Zbrozek  and  Jones  (Ref.  15).  In  Reference  16,  a  method 
for  estimating  the  response  of  a  wing  to  non-stationary  buffet  loads  was  developed  and  applied  to  a 
number  of  examples  involving  different  aircraft  maneuvers.  The  analysis  was  based  on  a  theory  of  the 
non-stationary  response  of  linear  dynamic  systems  first  given  by  Caughey  and  Stumpf  (Ref.  17),  and 
later  extended  by  Bamoski  and  Maurer  (Ref.  18),  and  Holman  and  Hart  (Ref.  19).  Various  forms  of 
the  power  spectral  density  for  the  random  part  of  the  input  load  were  investigated,  and  expressions 
were  derived  for  the  mean  square  response  of  the  wing  displacement.  In  Reference  20,  it  was  shown 
that  if  the  loading  on  the  wing  was  to  compose  of  fluctuating  forces  from  various  uncorrelated  sources, 
a  series  representation  could  be  used  to  represent  fairly  complex  buffet  load  spectra.  The  total  response 
could  be  obtained  from  a  summation  of  the  individual  responses  to  each  term  in  the  series.  The  analysis 
was  further  extended  in  Reference  21  to  include  multi-modal  response  of  a  wing  to  non-stationary 
buffet  loads  for  those  cases  when  the  undamped  natural  frequencies  of  individual  modes  were  close 
together  and  the  modes  of  vibration  were  statistically  dependent. 

This  report  describes  a  method  for  predicting  the  mean  square  displacement  and  acceleration 
response  of  a  wing  to  non-stationary  buffet  loads  by  extending  the  analyses  of  References  16,  20  and 
21 .  The  time  history  of  the  applied  load  is  segmented  into  a  number  of  time  intervals.  In  each  time 
segment,  the  non-stationary  load  is  represented  by  the  product  of  a  deterministic  shaping  function  and 
a  statistically  stationary  random  function.  The  total  response  to  a  given  maneuver  is  obtained  by  sum¬ 
mation  of  the  responses  of  the  wing  from  each  time  segment.  Unlike  the  approach  used  in  References 
16,  20  and  21,  which  assume  a  known  integrated  buffet  load  over  the  wing,  the  present  study  gives  an 


approximate  modelling  of  the  load  on  the  wing.  This  is  based  on  a  modification  of  a  method  proposed 
by  Schweiker  and  Davies  (Ref.  22)  for  the  study  of  the  response  of  shells  to  aerodynamic  noise.  In 
essence,  the  method  divides  the  wing  into  panels  or  elements,  and  the  load  is  computed  from  measured 
or  estimated  pressure  fluctuations  at  the  centre  of  each  panel.  Following  the  approach  used  in  Refer¬ 
ence  20,  a  series  representation,  with  terms  of  the  correlated  noise  type,  is  used  to  represent  the  complex 
buffet  pressure  power  spectral  densities.  Some  computed  response  results  are  presented  for  the  F-4E 
wing  as  an  example,  using  experimental  fluctuating  pressure  data  obtained  by  Mullans  and  Lemley 
(Ref.  9). 

2.0  ANALYSIS 

2.1  Dynamic  Aeroelastic  Equations 

Consider  a  Cartesian  co-ordinate  system  x,  y  and  z  fixed  on  the  wing  as  illustrated  in  Figure  1. 
The  displacement  of  the  wing  can  be  expressed  in  terms  of  a  set  of  normal  co-ordinates  Zj(t)  as: 

z(x,y,t)  =  2  0j(x,y)Zj(t)  (1) 

i  =  1 

where  0j(x,y)  is  the  mode  shape  function  of  the  ith  mode,  and  I  is  the  number  of  modes  required  to 
adequately  represent  z(x,y,t)  in  the  form  of  a  series.  The  dynamic  aeroelastic  equations  governing  the 
response  of  z((t)  to  an  input  load  Lj(t)  is  given  in  generalized  co-ordinates  as: 

Mjz’j  +  Cjij  +  KjZi  =  Lj(t)  (2) 

for  the  ith  mode.  The  dots  denote  differentiation  with  respect  to  time,  M,,  C,  and  K,  are  the  generalized 
mass,  damping  coefficient  and  stiffness  of  the  ith  mode  respectively.  In  the  above  equation,  the  assump¬ 
tion  of  light  damping  is  made,  thus,  cross  damping  terms  do  not  appear. 


Define  the  undamped  natural  frequency  coin  as: 


K, 

Mi 


and  the  damping  ratio  as: 


ci 

=  ITm^T 

Using  the  above  two  expressions,  Equation  (2)  can  be  rewritten  as  follows: 


+  2fiWinii  + 


Mi 


Lj(t) 


The  generalized  force  can  be  expressed  in  terms  of  the  fluctuating  pressure  p(x,y,t)  as: 


(3) 


(4) 


(5) 


Lj(t)  =  //  0j(x,y)p(x,y,t)dxdy  (6) 

wing 

planform 


where  the  integration  is  taken  over  the  wing  surface. 


2.2  Response  to  Non-Stationary  Input 

From  Equation  (1),  the  mean  square  response  of  the  displacement  of  the  wing  can  be  written 
in  the  following  form: 


E[z2(x,y,t)]  =  2  2  0i(x,y)0j(x,y)E[zi(t)zj(t)]  (7) 

i  =  l  j  =  1 

■y  -y 

For  i  =  j ,<j>.  (x,y)E[z.  (t) ]  gives  the  mean  square  response  of  the  ith  mode  at  the  point  (x,y)  on  the  wing. 

If  i  ¥=  j,  0i(x,y)0J(x,y)E[zj(t)Zj(t)]  is  the  contribution  to  the  total  response  due  to  correlation  between 
the  ith  and  jth  modes. 

For  a  given  buffet  maneuver,  the  mean  square  loading  on  the  wing  versus  time  may  be  re¬ 
presented  schematically  as  in  Figure  2.  The  points  A  and  B  correspond  to  the  onset  and  exit  of  buffeting 
respectively.  The  time  duration  the  aircraft  spends  in  the  buffet  regime  is  denoted  by  TB .  The  load  is 
non-stationary  and  time  varying  properties,  such  as  the  mean  square  of  the  response,  can  only  be  deter¬ 
mined  by  instantaneous  averaging  over  an  ensemble.  The  time  segmentation  technique  used  in  References 
16,  20  and  21  divides  the  time  TB  into  a  number  of  intervals  or  segments.  Within  each  segment,  the 
load  is  represented  by  the  product  of  a  deterministic  time  function  and  a  stationary  random  function. 
The  random  function  may  not  necessarily  be  the  same  in  each  segment.  Using  the  subscript  ‘r’  to  denote 
the  rth  segment,  the  generalized  force  for  the  ith  mode  at  time  t  can  be  written  as: 

Lj(t)  =  2  7r(t)fj  (t)  (8) 

r  =  l  1 

where  f,  (t)  is  a  stationary  random  function,  and  yr  (t)  is  a  deterministic  function  of  time  which  may  be 
considered  as  a  shaping  function  written  in  the  following  form: 


7r (t)  =  er0r(t) 

(9) 

where 

0,(t)  =  u(t-tr_, )-  u(t-tr) 

(10) 

Here  u(t)  is  a  unit  step  function,  er  is  a  constant  at  each  time  segment  and  represents  the  intensity  of 
the  generalized  input  load.  Figure  3  illustrates  the  variation  of  e  with  time  and  shows  an  approximate 
description  of  a  continuously  varying  load  using  the  time  segmentation  method.  However,  such  a 
representation  can  be  exact  in  practice,  for  example,  in  wind  tunnel  testing  where  the  model  angle  of 
incidence  advances  in  steps  with  time  in  a  prescribed  manner  so  that  the  buffet  load  can  be  represented 
by  Equation  (8)  without  any  approximation. 

Introduce  a  frequency  response  function  Hj(cj)  as 

Hj(cu)  =  — - - - - -  (11) 

Mj  J^jj2  -  to2  +  i2fjU>coinJ 

The  correlation  between  Lj(tj )  and  Lj(t2 )  can  be  written  as: 


n 

2 


r  =  1 


2  7f(ti  )7»(t2)^fi  fj  (^1 
>  =  I  *r  JI 


(12) 


-4- 

Since  f  (t, )  and  fj  (t2 )  are  assumed  to  be  stationary,  then 


f.  (t(  ,t2 )  -  Rf  f 
Js  r 


/  Sf.  f.  (cj)eiu>Tdto 

-oo  *r  h 


(13) 


where  t  -  t,  -  t2.  Using  Equations  (12)  and  (13)  and  the  generalized  Wiener-Khinchine  relation  (Ref.  23), 
SL .L.(coj  ,oj 2 )  can  be  expressed  as 


SLiLj(CJl  >w2 ) 


/  Sf.  f.  (co)Aj  (oo  -  co, ) A*  (oo  -  oo2)doo 
-oo  *r  1«  r  s 


(14) 


where  Ai(oo-oo,)  =  —  /  7r(t,  )ei(aj_“  l  )li  dt,  (15) 

r  2ir  _oo 

and  A*  (oo- oo2)  =  — -  /  7s(t2)e-i(“-“2>‘2dt2  (16) 

The  correlation  of  z,(t, )  and  Zj(t2 )  can  be  expressed  in  terms  of  the  power  spectral  density 
of  the  applied  load  SL  L,(oo,  ,oo2)  as  follows: 


RziZj(t,,t2)  =  f  J  SL.L.(oo1,oo2)Hi(oo1)H*(oo2)ei(wi,i-u’2,2>dco1doo2  (17) 


where  Hj(oo, )  is  given  by  Equation  (11 )  and  H*  (oo2 )  is  the  complex  conjugate  of  Hj(oo2  ).  Substituting 
Equation  (14)  into  the  above  equation  results  in  the  following: 


^ZiZ^i ’^2 )  ^  ^  /  Sf.  f.  (co)I,  (t,  ,oo)I;  (t2 ,oo)dco 

1  ‘  r  »  1  s  =  1  ’r  Js  r  JS 


(18) 


where  I(t,oo)  is  the  time  dependent  frequency  response  function  defined  as: 


I,r(t,,co)=  /  A,'(co  -  (o,)Hj(oo,  )e  * ' •  doo, 


(19) 


and 


Ij*s(t2,oo)  =  /  A*(co-co2)H*(co2)e-iw2‘2dco2 


(20) 


Assuming  f,(t)  and  fj  (t)  to  be  statistically  independent  between  any  two  time  segments, 
Equation  (18)  can  be  simplified  to  the  following: 


Rz ,z j^l  >^2 )  =  E(zi(ti  )zj(t2 )]  =  2  E[z,  (t,  )Zj  (t2 )] 

1  r  =  l  1  r 


(21) 
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where 


E(Zj  (t,)z,  (t2)]  =  /  &  f  (<o)I;  (t,  ,<o)I*  (t2  ,<o)dco 

r  r  _  c»  ‘r  Jr  *  r 


(22) 


Ij  (t,  ,co)  can  be  obtained  from  Equations  (15)  and  (19)  by  contour  integration  (Ref.  16). 
After  some  algebra,  the  final  form  is: 


lir(t,,w)  =  erHi(w)  *  |u(t,-tr_,)  [ei<J,i  -  eiwtr-i  (^(t,-^,)  +  ^(tj-t^j  )j] 

-u(t,-tr)  ^eUJti-e“otr  (^iftj-t,..  j )  +  i^-  <Mt,-tr))]J 


(23) 


where 


wid  =/l  - 


is  the  damped  natural  frequency  of  the  ith  mode,  and 


(coscV  +  fi“  sin  ^v) 


i//j(t)  =  e  Icosojj  t  +  f 


(24) 


^j(t)  =  e  sin  coj(1t 


(25) 


Similar  expressions  can  be  derived  for  IJr(t2  ,co)  and  hence  Ij*  (t2  ,co).  If  an  analytic  expression  for  the 
power  spectral  density  Sf.  f.  (co)  is  specified,  Equation  (22)  can  be  integrated  by  residue  calculus. 

'r  h 

Ij  (t,  ,co)Ij*  (t2  ,oj)  can  be  determined  from  Equation  (23)  and  it  is  given  in  Appendix  A. 


2.3  Approximate  Representation  of  Buffet  Loads  on  Wing 

In  the  previous  section,  it  is  seen  that  the  response  of  a  wing  to  non-stationary  buffet  loads 
can  be  evaluated  once  the  power  spectral  density  Sf.  f.  (to)  in  each  time  segment  is  known.  Let  pr(x,y,t) 

'r  h 

be  the  fluctuating  pressure  on  the  wing  at  the  rth  time  segment.  Here  pr(x,y,t)  is  taken  to  be  statistically 
stationary  and  is  related  to  the  function  fif  (t)  in  Equation  (8)  by  the  following  expression: 


fjr(t)  =  -  //  0j(x,y)pr(x,y,t)dxdy 

er 

wing 

planform 

The  power  spectral  density  can  be  written  as: 


(26) 


Sf.  f.  (oo)  =  —r  JUS  0i(Xi.yi)0j(X2-y2)Spr(x1*yi’X2’y2’W)dxldyldx2dy2 

'r  Jr  ’ 


(27) 


-6- 


where  Sp^Xj  ,y,  ,x2  ,y2  ,co)  is  the  cross  power  spectral  density  between  points  (x,  ,y, )  and  (x2  ,y2 ) 

(Fig.  1)  in  the  rth  time  segment. 

Following  a  method  used  by  Schweiker  and  Davis  (Ref.  22),  and  Davis  (Ref.  24)  for  the  study 
of  the  response  of  shells  to  aerodynamic  noise,  the  wing  is  divided  into  a  number  of  panels  (Fig.  4). 
Approximation  to  the  integral  in  Equation  (27)  is  made  by  assuming  the  pressure  at  the  centre  of  a 
panel  to  be  representative  of  the  pressure  field  at  all  points  within  that  panel.  That  is,  the  pressure 
spectrum  at  any  point  is  given  by  that  at  the  centre  of  the  panel.  The  pressure  cross  spectrum  between 
two  points  on  the  same  panel  is  taken  to  be  equal  to  the  pressure  spectrum  at  the  panel  centre,  while 
if  they  lie  on  different  panels,  it  has  the  same  value  as  the  cross-spectrum  between  the  centres  of  the 
two  panels.  Equation  (27)  can  be  simplified  to  the  following: 


K 

2 

k  =  1 


K 

2 

8  =  1 


f  f  Sp  (x,  ,y1x2,y2,oj)^i(x1  ,y, )  ^(xj ,y2 )dx,  dy,  dx2dy2 
Ak  Ae 

(28) 


where  K  is  the  number  of  panels  the  wing  is  divided  into,  Ak  and  Av  are  the  areas  of  the  kth  and  Cth 
panel  respectively.  Within  the  panel  of  area  Ak ,  an  average  mode  shape  for  any  mode  ‘i’  is  taken  as: 


-k  1 

4>i  =  —  /  0j(x,y)dAk  (29) 

Ak  Ak 

If  points  (x,  ,y, )  and  (x2  ,y2 )  fall  on  the  same  panel,  say  the  kth  panel,  S  (x,  .y,  ,x2  ,y2  ,co)  is  written 
as  SPr  (oj).  On  the  other  hand,  if  points  (x,  ,y, )  and  (x2  ,y2 )  lie  on  different  panels,  for  example,  the 
kth  and  Cth  panels,  then  Sp(x,  ,y,  ,x2,y2,w)  is  represented  by 


Spr(x,  ,yi  ,x2,y2,a>)  =  (<o)e 


(30) 


where  5k  is  the  spatial  decay  coefficient  measured  at  the  point  (xk  ,yk )  and  rks  is  the  distance  between 
the  centres  of  the  kth  of  £th  panels  (Fig.  4). 

Using  Equations  (29)  and  (30),  Equation  (28)  can  be  rewritten  as  follows: 


1  K 

Sfi  fj  2 

‘r  Jr  6f  k  = 


where 

Xk  =  0k  Ak  •  jtfk  Ak  + 


Using  Equations  (31)  and  (32)  to  determine  Sf.  f.  (cu),  and  upon  substitution  into  Equation  (22) 

*r  Jr 

and  carrying  out  the  integration,  the  expected  value  E[Zj(t,  )zj(t2 )]  can  be  evaluated  from  Equation  (21) 
using  Equation  (22). 


Skr(oo)Xk 


k 


<l>*  Aee 


-6 


krke 


(31) 


(32) 


m 


i 


For  a  given  wing,  0k  and  0k  can  be  calculated  along  with  Ak  and  rkc  once  the  dimensions 
of  the  panels  have  been  decided.  Sk  (cj)  and  5.  are  usually  determined  experimentally  or  are  estimated 

P  T  K 

based  on  available  power  spectral  density  and  decay  coefficient  data  on  similar  wings.  There  are  certain 
forms  of  the  pressure  spectrum  that  enable  the  integral  in  Equation  (22)  to  be  evaluated  analytically. 
The  pressure  at  any  point  on  the  wing  can  be  considered  to  be  made  up  of  contributions  from  a  number 
of  sources.  These  sources  may  arise  from  flow  separations,  wing  leading  edge  and  tip  vortices  and  other 
disturbances  characteristic  of  that  particular  wing  geometry.  Assuming  that  the  pressure  fluctuations  of 
these  sources  are  uncorrelated,  Sk  (u)  can  be  written  in  series  form  as: 


M 

Sk  (co)  =  2  SkmM 

r  m  =  1 


(33) 


where  M  is  the  number  of  terms  that  give  a  good  fit  to  the  experimental  data.  Reference  21  shows 
that  quite  complex  power  spectral  density  curves  can  be  derived  from  the  series  representation  of 
Equation  (33)  using  the  correlated  noise  and  exponential  (commonly  used  in  isotropic  turbulence) 
froms  of  the  power  spectral  density. 

2  4  Mean  Square  Response  of  Wing  Displacement  for  Given  Input  Power  Spectra 

From  Equation  (21),  the  expected  value  Efz^tj  )z.(t2)]  can  be  rewritten  as: 


E  [Zi(tj  )Zj(t2 )]  =  E  [ziq  (t,  )Zj^  (t2  )J  +  ‘V  E^t,)^)]  (34) 


where  tq_ ,  <  tj  <  tq  and  tq_  j  <  t2  <  tq .  In  other  words,  if  t,  and  t2  lies  in  the  qth  time  segment, 
the  expected  value  is  given  by  the  sum  of  E[Zj  (tj  )Zj  (t2 ))  at  the  qth  segment  and  contributions  from 

the  r  =  1  to  q  -  1  segments.  The  significance  of  writing  E[Zj(t,  )zj(t2 )]  in  the  form  of  Equation  (34) 
will  become  clear  later  on  when  it  is  shown  that  the  two  terms  on  the  right  hand  side  of  the  equation 
are  obtained  from  different  expressions. 

If  Skm(w)  takes  on  the  spectrum  of  a  white  noise,  then  Equation  (33)  can  be  written  as: 

P  r 


Sk  («)  =  Sk  =  constant 

pr'  '  r 


(35) 


Substituting  Equations  (31)  and  (35)  into  (22)  yields 


E [»,,(*. >V‘2>]  'if  kr, x‘£[vi, 


(36) 


fi  [z,r(t,  )Zj  (t2)]  =  f  SkI.r(t,  ,w)I*  (t2,w)dw 


where 


(37) 
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Equation  (37)  can  be  evaluated  for  tr_ ,  <  t,  <  tr ,  tr_  ]  <  t2  <  tr  using  Equation  (Al)  of  Appendix  A, 
while  Equation  (A2)  is  used  for  t,  >  tr,  t2  3*  tr.  The  poles  in  the  integrand  are  given  by  the  term 
HiMH'  (cu),  and  there  are  four  of  them  in  the  co-  plane,  namely,  p,  =  coid  +  ifj0Jin  .  P2  = 
p-,  =  co:  -  if; co:  and  p4  =  -co; ,  -  ifjCJ;  .  Using  residue  calculus,  the  integrations  can  be  performed 
by  choosing  the  proper  integrated  paths.  The  resulting  expressions  for  E  [z,  (t,  )z]  (t2 )]  are  given  by 
Equations  (Bl)  and  (B2)  in  Appendix  B. 

The  mean  square  response  of  the  wing  displacement  is  obtained  from  Equation  (7)  using 
Equation  (34)  after  setting  t(  =  t2  =  t.  If  Equations  (36)  and  (37)  are  substituted  into  Equation  (34) 
to  obtain  E[Zj(t]  )zj(t2 )],  the  resulting  response  is  for  a  white  noise  pressure  spectrum  in  each  of  the 
time  segments. 

To  handle  more  complex  buffet  load  spectra,  the  term  Sk  m(co)  in  Equation  (33)  can  be 
written  in  the  following  form: 


S£m(o;)  = 


Skma2 


km 


(*km  .km  ,  \ 

0i\  +  02  +  co2) 


[a2  k  m  +  (co  +0[cm)2J  [a2 


,  km  ,  2l 

2  +(co-0rkm)J 


(38) 


which  is  the  correlated  noise  power  spectral  density.  Here  S,  a  and  0  are  constants,  and  the  subscript  ‘r’ 
denotes  the  rth  time  segment,  the  superscripts  ‘k’  and  ‘m’  denote  the  kth  panel  and  the  mth  mode 
respectively.  Note  that  if  0km  is  set  equal  to  zero,  the  resulting  expression  is  the  same  as  that  used  in 
Reference  16  for  power  spectrum  of  the  form  used  for  isotropic  turbulence. 

Using  Equations  (31)  and  (33),  Equation  (22)  becomes 


where 


E  [zj  ft,  )Zj  (t2)]  =■— 
Ekn,[  z.  (ti  )zj  (t2)]  = 


K  M  .  ^  r  1 
2  2  xkEkm  I  1  (t,  )z.  (t2 ) 

k  =  1  m  =  1  L  r  1  Jr  1  -J 

/  Skm  (co)Ij  (t,  ,co)I*(t2  ,co)dco 

00  r  r 


(39) 


(40) 


Substituting  Equations  (Al)  and  (A2)  into  the  above  expression,  and  using  Equation  (38)  for 
S*5™  (co),  £km  [Zj(t|  )zj(t2)]  for  tr_ !  <t,  <tr,tr_!  <t2  <  tr  is  given  by  Equation  (Cl).  When  t,  3*  tr, 
t2  3s  tr,  the  expected  value  is  given  in  Equation  (C2).  The  mean  square  value  of  the  wing  response  is 
again  given  by  Equations  (7)  and  (34)  with  t,  =  t2  =  t  and  using  Equations  (39)  and  (40)  to  obtain 
E[zif(t)zjf(t)]. 


2.5  Mean  Square  Response  of  Wing  Acceleration  for  Given  Input  Power  Spectra 

An  expression  similar  to  that  of  Equation  (7)  can  be  written  for  the  acceleration  by  simply 
replacing  z(x,y,t)  with  the  second  time  derivative  z  (x,y,t),  that  is, 


A 


E[z2(x,y,t)]  =  2  2  0i(x,y)0j(x,y)E[zi(t)zj(t)] 

i  =  i  j  =  i 


E[Zj(t)z.(t)]  is  obtained  from  the  following  equation: 

q  -  l 

EQ^t^t)]  =  E[k\  (t)iij  (t)]  +  2  ^  E[zif(t)zjr(t)]  (42 

where  tq_  t  <  t  <  tq .  If  the  pressure  fluctuations  in  each  time  segment  have  spectra  of  the  white  noise 
type,  then: 


1  J? 


E[zj  (t)Zj  (t)]  =  — r-  2  XkE[zi  (t)Zj 

t  Ji  J  ef  k  =  i  r 


while  for  correlated  noise  power  spectral  density. 


1  K  M 


E[z.  (t)Zj  (t)]  =  — r  2  2  XkEkm[zj  (t)z.  (t)J 

r  Jr  €•  k  =  1  m  =  I  r  Jr 


where 


e[z;  (t)z\  (t)]  =  — - — —  £[z.  (t,)zj  (t2)] 
1  r  dt2at2 


t,=t2  =  t 


skm[zt  (t)z  (t)] - Ekm[zj  (t  )z  (t2)] 

1  Jt  atfatf  r  Jr 


1 1  =  *2“ ' 


Using  the  equations  given  in  Appendices  B  and  C,  Equations  (45)  and  (46)  can  be  evaluated 
by  carrying  out  the  differentiations  and  setting  tj  =  t2  =  t.  The  resulting  expressions  are  substituted 
into  Equations  (43)  and  (44).  The  mean  square  response  of  the  acceleration  of  the  wing  is  then  de¬ 
termined  from  Equations  (41)  and  (42).  The  expressions  for  £[Zjr(t)Zj  (t)]  are  given  in  Appendix  D, 

while  those  for  Ekm  [zjf (t)z^(t)]  are  given  in  Appendix  E. 


2.6  Expressions  for  Displacement  and  Acceleration  Response  in  Non-Dimensional  Forms 

It  is  convenient  for  computation  purposes  to  express  the  equations  derived  in  Sections  2.4 

and  2.5  in  non-dimensional  forms.  Consider  first  the  case  where  the  pressure  power  spectral  densities 

are  of  the  form  for  white  noise.  Panel  1  (k  =  1)  is  taken  as  the  reference  panel  and  the  subscript  ‘rQ  ’  is 

used  to  present  the  conditions  at  e(t)  =  1.  Using  the  time  segmentation  technique  (Fig.  3),  AC„ 

Prms 

and  the  power  spectral  density  curves  given  in  Reference  9  are  assumed  to  be  the  values  at  the  time 
segment  when  the  load  intensity  is  at  its  maximum,  that  is,  at  er  =  1 .  The  relation  between  Sk  at  the 
rth  time  segment  and  its  value  at  er  =  1  is 

Sk  =  er2Sk  (47 

f  i  i 
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Define  the  following  non-dimensional  quantities: 


In  terms  of  Sk  ,  Sk  can  be  written  as: 


Sk 

r 


(48) 


k  - 


(49) 


<t>* 


k 


*J 


(50) 


A,  =  — 


(51) 


f. 


fi 

FT 


(52) 


M. 


M. 

I 

mT 


(53) 


co.  — 

‘n  CO, 


(54) 


and  S  k  can  be  expressed  as: 


sf  -  e]  sf 


sf 


fK. 


WCJ. 


(55) 


(56) 


From  Equation  (32),  let 


Xk  -  Xk(#})2  Aj  2 


(57) 


where  Xk  =  0kAk{0kAk  +  2  0?Aee6krkcj 

Substituting  into  Equation  (36)  and  using  either  Equation  (Bl)  or  (B2)  for  ElZj  ftj  )Zjr(t2 )]  yields  the 
following  expression  for  E[Zj  (t)Zj  (t)]: 

tf'F  i  k  — .  _ 

E[zj  (t)z ,  (t)]  =  e2(0>)2A2KSj  - -  2  XkSk  (58) 

'  J'  r  af-MjMjwf  K  k  =  l 


In  deriving  this  expression,  t(  and  t2  have  been  set  equal  to  t.  Let  EfZj  (t)zj(t)]  be  the  non-dimensional 
form  of  E[zit(t)Zj  (t) ]  defined  as: 


E[z.  (t)z  (t)]  = 

t  Jr 


2f.  M2gu* 

n 

ttK(0})2A2sJ 


-E[Zj  (t)z.  (t)] 
r  Jr 


Then  E[Z(  (t)Zj  (t)]  can  be  expressed  as: 


e2  'F„ 


E[Zj  (t)z,  (t)]  =  _  ■  r  X 

r  K  k  -  i 


'  '*  X  x‘sf 


(59) 


(60) 


where  the  expressions  for  >Fa  for  tr_ ,  <  t  <  tr ,  and  t  >  tr  are  given  in  Appendix  B.  A  similar  equation 
can  be  derived  from  Equation  (39)  for  correlated  noise  power  spectral  densities.  Since  it  is  assumed 
that  the  shapes  of  the  power  spectral  density  curves  at  different  time  segments  are  similar  and  differ 
in  scale  only,  a*  m  and  /3k  m  in  Equation  ( 38 )  can  be  written  in  terms  of  the  values  at  er  =  1  as  m 

and  0km  respectively. 

ro 


Define 


Sk  (0) 

F  r 


(61) 


and  similar  to  Equation  (47),  let 
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Using  Equation  (39)  and  Equation  (Cl)  for  tr_  [  <  t  <  tt,  or  Equation  (C2)  for  t  >  t,,  Elzjt(t)z^(t)l 
can  be  written  as: 


E[Zj  (t)z.  (t)]  = 

I  JI 


K  M  . 

--2  2  XkSkm'I' 


2fjMjMjcj|3  k  =  l  m  =  l  rc 


(63) 


Using  Equation  (57),  define 


'AC_ 


Xk 


r  m» 


i  f  M 

Acp  /tt  2  §kmakm 

F  rms/  r  r 

m  =  l  0  ° 


(64) 


and 


M 


s0  -  7T(0|)2a2s'  (0)  2  S  ’  m  Q,  1  m 

rn  m  =  1  0  0 


(65) 


Similar  to  Equation  (59),  define  a  non-dimensional  form  for  ElZj  (t)Zj  (t)]  as  follows: 


E[Zj  (t)Z:  (t)]  = 
r  Jr 


2?jM2co3 


ttKS„ 


E[z4  (t)Zj  (t)] 
r  Jr 


(66) 


Combining  Equations  (63)  to  (66)  yields. 


E[z;  (t)z,  (t)]  =  z~~ZT~  ~  2  2  XkSkm^b 

r  Jf  fjM.M.w?  K  k  •  •  • 

1  3  n 


m  =  1 


(67) 


The  acceleration  response  for  white  noise  form  of  pressure  power  spectral  density  can  be 

obtained  from  Equation  (42).  If  the  non-dimensional  form  for  E[z.  (t)z,  (t)]  is  written  as: 

‘r  Jr 


-  ..  ..  *.“?  ..  .. 

E(Zj  (t|z  (tt)  - - _•  :  —  Elzi  W*.  <l)l 

r  Jr  w,  wK(0  j  )2  A2Sj  J' 


(68) 


then  using  Equation  (Dl)  for  tr_  j  <  t  <  tr,  or  Equation  (D2)  for  t  >  tr,  it  can  be  expressed  as: 


E[Zj  (t )*,  (t)]  =  er2^  2  XkSk 

r  1  CK  k  =  l 


(69) 
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Corresponding  to  Equation  (44)  for  correlated  noise  power  spectral  density,  if  E  [z,  (t)zj  (t)]  is  defined 
as:  r  ' 


-  ..  ..  2fiM?  ..  .. 

E [z  (t)z.  (t)]  =  - ;E[Zj  (t)zs  (t)j 

‘r  Jr  a1 1  ttKS  J' 

ro  0 


(70) 


then  using  Equations  (El)  or  (E2)  gives  the  following: 


E[Zj  (t)Zj  (t)]  =  — 

1 


1  K 
—  2 
K  k  =  1 


M 

2 

m  =  1 


a' 


km 


xkskm*d 


a 


li 


(71) 


where  S0  is  given  by  Equation  (65).  With  expressions  E  [zj  (t)Zj  (t)]  and  E  [z,  (t)Zjf(t)|  given  by  Equa¬ 
tions  (60)  and  (68),  or  Equations  (67)  and  (71)  for  white  noise  or  correlated  noise  form  of  power  spectral 
densities  of  the  pressure  fluctuations  respectively,  the  total  displacement  and  acceleration  response  at 
any  point  on  the  wing  can  be  obtained  from  Equations  (7)  and  (41).  Non-dimensional  forms  of  the  mean 
square  response  E  [z2(x,y,t)]  or  E  [z2(x,y,t)l  can  be  obtained  following  the  definitions  of  E  (Zi  (t)Zj  (t)] 
and  E  [zif  (t)zJr(t)]  given  in  Equations  (59)  and  (66),  or  Equations  (68)  and  (70)  respectively. 


3.0  AN  ILLUSTRATIVE  EXAMPLE 

The  steps  involved  in  computing  the  response  of  the  wing  are  fairly  straightforward  once  the 
necessary  aeroelastic  data  are  given.  As  an  example,  the  response  of  a  F-4E  wing  is  studied  using  the 
data  supplied  by  Mullans  and  Lemley  (Ref.  9).  Figure  5  is  taken  from  Reference  9  and  shows  the  choice 
of  panel  dimensions  on  the  wing.  A  total  of  eighteen  panels  (K  =  18)  are  used  in  the  computations,  and 
the  panel  centres  are  marked  numerically  from  1  to  18.  The  distances  between  the  centres  of  any  two 
panels  rke  can  be  determined  from  the  geometry  of  the  wing  planform. 

Only  the  first  ten  symmetrical  modes  are  considered  in  this  example.  The  natural  frequencies 
:  ,  generalized  mass  M=,  average  mode  shape  0  k  and  damping  ratio  for  dynamic  pressure  q  =  470  psf 

are  tubulated  in  Mullans  and  Lemley’s  report  (Ref.  9).  They  also  obtained  buffet  data  from  wind  tunnel 
tests  on  a  10%  scale  rigid  three  dimensional  model.  For  a  clean  wing  with  zero  leading  edge  and  trailing 
edge  deflections,  the  rms  pressure  coefficients  ACL  on  the  wing  are  shown  in  Figure  6.  In  this  par- 

'  rmi 

ticular  example,  AC_  varies  practically  linearly  along  constant  chord  and  span  directions. 
v  rm  s 

The  power  spectral  densities  of  the  pressure  fluctuations  at  various  points  on  the  wing  for 
Mach  number  0.7  and  angle  of  incidence  12°  are  also  given  in  Reference  9.  The  transducers  are  not 
located  exactly  at  the  panel  centres  and  there  may  be  more  than  one  transducer  per  panel,  or,  in  some 
cases,  none  at  all.  Figure  7  indicates  the  transducers  (numbers  correspond  to  these  given  in  Ref.  9) 
whose  measured  outputs  are  used  to  represent  the  fluctuating  pressures  on  the  panels  where  they  are 
located. 


The  value  of  the  spatial  decay  coefficient  8k  given  in  Equation  (30)  can  be  found  from 
cross-correlation  of  the  pressure  fluctuations  at  two  points  on  the  wing.  Its  value  is  taken  to  be  con¬ 
stant  between  any  two  points  on  different  panels  and  equal  to  0.25.  The  variation  of  the  intensity  of 
the  input  load  with  time  is  assumed  to  be  sinusoidal,  that  is, 

t 

e(t)  =  sin  it  — 

tb 
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where  T„  is  the  time  the  aircraft  spends  in  the  buffet  regime  (Fig.  3).  The  values  of  ACn  given  in 

v  rmi 

Figure  6  are  taken  to  correspond  to  e(t)  =  1,  or  when  t  =  TB  /2.  Furthermore,  it  is  assumed  that  the 
normalized  shapes  of  the  pressure  power  spectral  density  curves  with  respect  to  its  value  at  zero  fre¬ 
quency  at  different  panels  are  the  same  for  all  values  of  e(t).  In  other  words,  varying  e(t)  changes  the 
magnitude  of  power  spectra  only,  but  the  normalized  values  are  the  same  for  all  time  segments.  This 
requires  the  flow  field  on  the  wing  to  be  similar  at  all  times  during  the  buffet  maneuver  with  changes 
only  in  the  scale  of  the  fluctuating  pressure  force.  This  should  be  viewed  as  a  rough  approximation 
only.  A  more  complete  representation  requires  measurements  of  ACn  and  the  pressure  power 

*  r  ms 

spectral  densities  on  all  panels  at  flight  conditions  which  correspond  to  the  value  of  er  under  consid¬ 
eration. 


4.0  RESULTS  AND  DISCUSSIONS 


In  the  formulation  given  in  the  previous  sections,  the  Mach  number,  angle  of  incidence  and 
dynamic  pressure  do  not  appear  explicitly  in  the  equations  defining  the  response  of  the  wing.  The 
input  buffet  data  necessary  for  the  computations  are  the  AC„  and  the  power  spectral  densities 

of  the  pressure  fluctuations  at  each  of  the  18  panels  for  the  time  segment  er  =  1.  The  data  used  in  this 
example  at  ex  =  1  correspond  to  those  given  by  Mullans  and  Lemley  (Ref.  9)  for  q  =  470  psf,  M  =  0.7 
and  a  =  12°.  The  variation  in  the  buffet  load  with  time  arises  from  changes  in  flight  conditions,  for 
example,  changes  in  the  angle  of  incidence.  A  sinusoidal  variation  of  e(t)  is  chosen  for  convenience 
in  the  computations,  even  though  the  time  segmentation  technique  can  be  used  for  quite  general 
buffet  maneuvers. 


The  experimentally  determined  power  spectral  density  curves  are  expressed  in  series  form 
using  Equations  (33)  and  (38).  The  unknowns  Skm ,  a,"1  and  0jcm  and  the  number  of  terms  in  the  series 

0  o  o 

are  obtained  by  a  curve  fitting  procedure  using  Powell’s  minimization  algorithm  (Ref.  25).  Figure  8 

(«)  (k  =  18) 


shows  a  comparison  between  the  curve  fitted  non-dimensional  power  spectral  density  S  £ 


versus  frequency  cu  for  transducer  16  at  panel  18  (Fig.  7)  with  the  smoothed  experimental  data  taken 
from  Reference  9.  It  is  seen  that  only  three  terms  are  necessary  in  Equation  (33)  to  give  a  very  good 
fit  to  the  experimental  data.  The  individual  terms  are  plotted  in  Figure  9.  Curve  ‘a’  has  the  exponential 
form  of  power  spectral  density  which  is  used  commonly  for  isotropic  turbulence,  and  has  been  inves¬ 
tigated  in  detail  in  Reference  16.  Curves  ‘b’  and  ‘c’  have  power  spectra  of  the  correlated  noise  type 
which  has  been  studied  in  Reference  20.  On  examining  these  curves,  the  pressure  can  be  decomposed 
into  three  components:  a  fluctuating  part  due  to  flow  separations  with  power  spectral  density  represented 
approximately  by  the  exponential  form,  and  two  convected  vortices  or  disturbances  with  correlated 
noise  power  spectral  densities.  These  convected  vortices  can  be  identified  as  the  leading  edge  and  wing 
tip  vortices.  In  curve  fitting  the  power  spectral  densities  taken  from  Reference  9,  it  is  found  that  three 
terms  are  needed  at  most.  For  some  panels,  two  terms  are  sufficient. 


The  mean  square  displacement  response  E[z2(x,y,t)]  given  in  Equation  (7)  can  be  written 
in  non-dimensional  form  at  a  panel  centre  as: 


Ek[z2(t)] 


1 

2 

i  =  1 


EfZjOdZjft)] 


(73) 


E[z;(t)Zj(t)]  is  obtained  from  Equation  (34)  by  setting  t,  =  t2  =  t  and  using  non-dimensional  forms  of 
the  expected  values.  Consider  e,  to  be  constant,  and  without  loss  in  generality,  take  er  =  1.  This  is 
equivalent  to  the  case  of  the  wing  response  to  one  time  segment  only,  and  Equation  (67)  can  be  used 


in  place  of  Equation  (34)  to  evaluate  E[Zj(t)Zj(t)).  The  term  in  Equation  (Cl)  for  tT_ ,  <  t  <  tr. 


after  setting  t,  =  t2  =  t,  determine  the  response  to  step  modulated  buffet  loads,  while  for  t  >  tr, 
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substitution  of  Equation  (C2)  into  Equation  (67)  gives  the  decay  to  pulse  modulated  loads.  The  pulse 
duration  is  given  by  the  time  segment  Atr,  which  is  equal  to  tr  -  tr_ , . 

Figure  10  shows  the  response  at  the  centre  of  panel  18  for  er  =  1  using  the  values  of  ACp 
and  spectral  information  given  in  Reference  9.  The  superscript '  is  used  to  denote  that  time  is  non- 

2ir 

dimensionalized  with  respect  to  the  undamped  natural  period  of  the  first  mode  of  vibration  T.  =  —  . 

03, 

1  n 

The  mean  square  response  is  normalized  with  respect  to  E  k  [z2  (t')]max ,  which  is  the  maximum  value 
obtained  as  t'  -  tj_  j  -*  °°.  This  maximum  response  is  the  sum  of  the  responses  from  the  ten  individual 
symmetrical  modes  of  vibration  together  with  the  contributions  from  the  correlation  terms  between 
different  modes.  The  magnitude  of  the  correlation  terms  can  be  obtained  from  the  difference  between 
curves  ‘1’  and  ‘2’  which  amounts  to  approximately  8%  for  t'  -  t'_ ,  >  1.  As  expected,  the  first  bending 
mode  of  vibration  is  the  dominant  one,  but  it  only  accounts  for  55%  of  the  total  response.  Because 
the  damping  is  fairly  large  for  this  particular  wing,  the  response  reaches  its  maximum  in  approximately 
three  undamped  natural  periods  Tj  of  the  first  mode.  The  oscillations  in  the  response  curves  are  rela¬ 
tively  small  as  compared  to  those  of  Reference  16  for  small  values  of  damping. 

The  decay  of  the  mean  square  displacement  response  at  the  centre  of  panel  18  after  the 
exciting  force  has  been  withdrawn  at  time  t'  is  shown  in  Figures  11  to  13  for  three  values  of  the  pulse 
duration  (At'  =  0.5, 1  and  2).  The  shapes  of  curves  are  very  similar  for  the  three  values  of  At'  shown. 

In  less  than  three  undamped  natural  periods  T,  of  the  first  mode,  the  decay  curve  is  practically  zero. 
Comparison  of  curves  ‘1’  and  ‘2’  shows  the  correlation  terms  to  be  small,  especially  for  increasing  At' . 
The  oscillating  behaviour  of  the  decay  curves  can  be  observed  from  these  figures,  but  the  amplitudes 
are  much  smaller  than  those  shown  in  Reference  16  for  small  damping  ratios. 

Figure  14  shows  the  response  to  a  step  modulated  input  buffet  load  at  the  centre  of  panel  18 
assuming  white  noise  pressure  power  spectral  densities  for  all  the  panels.  In  this  case,  Equation  (60)  is 
used  in  place  of  Equation  (67).  Comparison  with  Figure  10  shows  the  response  curves  to  be  nearly  the 
same.  Similar  observations  are  also  made  when  the  decay  curves  given  in  Figures  15  to  17  are  compared 
to  Figures  11  to  13.  It  thus  appears  that  approximating  the  power  spectral  densities  of  the  fluctuating 
pressures  in  all  the  panels  by  that  of  a  white  noise,  the  shapes  of  the  response  and  decay  curves  when 
one  time  segment  is  considered  do  not  differ  much  from  those  computed  using  series  representation 
of  the  experimentally  determined  power  spectral  densities.  It  has  been  shown  in  Reference  16  that 

for  m  =  1  and  =  0,  the  response  and  decay  curves  are  similar  to  those  with  white  noise  power 
spectral  densities  if  c c  is  small  compared  to  akm .  For  3km  =£  0,  which  corresponds  to  the  case  where 
the  power  spectral  density  is  of  the  correlated  noise  type,  the  results  presented  in  Reference  20  show 
that  if  03 i  is  not  near  0km ,  the  response  and  decay  curves  for  co;  /akm  <  1  are  similar  to  those  obtained 

*n  *o  *n  *o 

in  Reference  16  using  the  white  noise  form  of  power  spectral  density.  The  series  representation  of  the 
experimental  data  shows  that  values  of  (m  >  2)  are  large  compared  to  the  first  few  natural  fre¬ 
quencies  of  the  modes  that  contribute  significantly  to  the  response  and  decay.  Also,  at  those  panels 

where  values  of  ACD  are  large,  03.  <  akm  for  the  lower  modes.  Hence,  the  use  of  white  noise 
nmi  ‘n  ro 

power  spectral  densities  for  the  pressure  fluctuations  is  a  fairly  good  approximation  for  this  particular 
example.  It  should  also  be  noted  that  the  computation  time  using  white  noise  representation  is  signifi¬ 
cantly  less  than  that  using  series  representation.  A  comparison  between  Equations  (60)  and  (67)  shows 

K 

that  4'.  has  to  be  computed  X  Mk  times,  where  Mk  is  the  number  of  terms  required  in  the  series 

k  =  l 

expression  for  Sk  (w)  in  the  kth  panel,  while  in  Equation  (60),  ♦  needs  to  be  evaluated  only  once. 
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Figure  18  shows  the  maximum  displacement  response  (t'  -  t'_ ,  -*  °°)  normalized  with  respect 
to  that  at  panel  18  for  all  the  panels.  The  numerals  1  to  18  designate  the  panel  centres,  and  their  loca¬ 
tions  are  given  by  y/b.  Here  b  is  the  semi-span,  and  y  is  the  distance  measured  from  the  wing  root  along 
a  constant  percentage  chord  line  (Fig.  5).  As  expected,  the  response  increases  in  the  outboard  direction, 
and  is  larger  for  panels  near  the  trailing  edge  than  those  close  to  the  leading  edge. 

For  an  input  buffet  load  with  sinusoidal  variation  of  e  with  time,  the  displacement  response 
versus  time  t'  is  shown  in  Figure  19.  The  results  are  for  panel  18  using  white  noise  pressure  power  spectral 
densities  for  all  the  panels.  The  duration  of  the  applied  load  Tg  is  forty  times  the  natural  period  T,  for 
the  first  vibration  mode.  Two  cases  are  given  for  At'  =  0.5  and  2,  which  correspond  to  eighty  and  twenty 
time  segments  respectively.  Similar  to  Reference  16,  it  has  been  found  that  changing  the  number  of  time 
segments  produces  practically  no  differences  in  the  computed  results  if  At|  <  0.5.  In  all  subsequent 
computations,  At'  is  taken  to  be  0.5. 

It  has  been  shown  in  Reference  16  that  for  rapid  maneuvers,  that  is,  Tg  is  small,  the  maximum 
response  lags  the  applied  load,  and  the  amplitude  is  lower  than  that  obtained  if  Tg  is  very  long.  Since  the 
damping  for  the  wing  in  this  example  is  fairly  large,  the  lag  is  found  to  be  negligibly  small.  The  displace¬ 
ment  amplitude  decrement,  which  is  defined  as 

Ek[z2(t')l  ,  T .  —  E  k  [z2  (t')] 

5  .  - - - 5 -  X 

Ek[z2(t')]max  for  „ 

is  plotted  versus  Tg  in  Figure  20  for  panel  18.  6  decreases  sharply  with  Tg  initially,  and  it  is  less  than 
1%  after  14  periods  of  vibration  in  the  first  mode. 

The  mean  square  acceleration  response  to  step  and  pulse  modulated  buffet  loads  has  been 
computed  using  an  expression  similar  to  Equation  (73)  with  z(t)  replaced  by  z(t).  Comparing  the  results 
suits  for  white  noise  and  those  using  series  representation  for  the  experimental  pressure  power  spectral 
densities  yields  similar  conclusions  as  with  the  displacement  response.  Only  the  white  noise  power 
spectral  densities  results  are  presented  in  this  report,  and  they  are  shown  in  Figures  21  and  22  for 
panel  18.  The  results  indicate  that  the  response  reaches  its  maximum  value  in  less  than  one  natural 
period  T,  of  the  first  mode.  The  contributions  from  the  correlation  terms  are  larger  than  those  for 
the  displacement  response,  reaching  a  value  of  nearly  20%  of  Ek  [z2(t')]max . 

The  decay  curves  in  Figure  22  are  shown  for  two  values  of  the  duration  of  the  time  segment 
(At'  =  0.5  and  1).  The  curves  reach  practically  zero  value  for  t'  -  t'  ~  1.  As  in  the  response  case,  the 
correlation  terms  between  modes  are  more  significant  than  those  found  for  the  displacement.  For  the 
two  values  of  At|  shown  in  the  figure,  the  curves  are  practically  identical.  The  mean  square  accelera¬ 
tion  response  with  sinsoidal  variation  of  e  with  time  is  shown  in  Figure  23  for  Tg  =  40,  At'  =  0.5.  For 
the  same  value  of  At' ,  Figure  19  shows  the  oscillations  in  the  displacement  response  curve  to  be  much 
smaller. 
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5.0  CONCLUSIONS 

A  method  for  predicting  the  mean  square  displacement  and  acceleration  response  of  a  wing 
to  non-stationary  buffet  loads  has  been  developed  and  applied  to  an  example  using  available  data  for 
the  F-4E  wing.  The  time  segmentation  technique,  which  is  an  approximation  to  a  continuously  varying 
load  with  time,  gives  results  which  are  practically  independent  of  the  duration  of  the  time  segment 
when  it  is  below  a  certain  minimum  value.  In  this  particular  example,  the  minimum  value  is  found  to 
be  about  half  the  undamped  natural  period  of  the  first  vibration  mode. 


The  present  method  of  modelling  the  buffet  load  on  the  wing  may  overestimate  the  response, 
since  the  pressure  cross-spectrum  between  two  points  on  the  same  panel  is  taken  to  be  equal  to  the 
pressure  spectrum  at  the  panel  centre.  Also,  for  two  points  lying  on  different  panels,  it  has  the  same 
value  as  the  cross-spectrum  between  the  centres  of  the  two  panels.  This  simplification  is  probably 
admissible  in  providing  approximate  response  data  in  the  early  stages  of  aircraft  design,  since  detailed 
correlation  measurements  are  extremely  time  consuming  and  expensive  to  carry  out. 

The  complex  experimental  buffet  pressure  spectra  can  be  expressed  in  the  form  of  a  series  by 
a  curve  fitting  procedure.  The  pressure  can  be  decomposed  into  components.  It  is  found  that  for  all 
panels,  the  first  term  in  the  series  representation  gives  the  exponential  form  of  power  spectral  density, 
which  can  be  considered  to  be  a  fair  approximation  for  separated  flows.  Other  terms  indicate  the  pre¬ 
sence  of  convected  vortices,  and  the  number  of  such  vortices  depends  on  the  location  of  the  panel.  In 
this  particular  example,  it  is  found  that  approximating  the  pressure  power  spectral  densities  at  each 
panel  by  that  of  a  white  noise  yields  results  which  are  fairly  close  to  those  obtained  using  the  experi¬ 
mentally  measured  values. 

A  significant  portion  of  the  wing  response  can  be  accounted  for  by  the  first  few  lower  vibra¬ 
tion  modes.  The  correlation  between  modes,  which  are  considered  to  be  statistically  dependent,  is 
investigated  and  found  to  contribute  less  than  10%  to  the  total  mean  square  displacement  response. 

For  the  acceleration  response,  the  correlation  terms  are  nearly  20%  of  the  total. 

A  sinusoidal  variation  of  load  intensity  with  time  is  used  as  an  example  to  demonstrate  the 
computation  of  the  response  for  the  particular  wing  under  investigation.  Since  the  damping  for  this 
wing  is  fairly  large,  the  lag  between  the  maximum  response  and  the  applied  load  is  negligibly  small. 

The  displacement  amplitude  decrement  is  practically  zero  after  a  few  undamped  natural  periods  of 
the  first  vibration  mode. 
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The  expressions  for  ^  and  are  given  by  Equations  (A3)  and  (A4). 
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=  Hlffl 

~  H27J2 
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n 

- I2V4 

akm  2  * 

r 

(E12) 

to. 
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H, ,  H2 , 1,  and  I2  are  given  in  Appendix  C,  and 


T7,  =  —  (1  -  4f2)  (E16) 

t?2  =  ?j(3  -  4$f)  (E17) 

Wid 

‘'i  =  zr  (1 "  4?/}  (E18) 

in 

v2  =  fj(3  -  4ff)  (E19) 
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